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Abstract 

We present observations of the embedded massive young stellar object (YSO) candidate 08576nr292, 
obtained with X-shooter and SINFONI on the ESO Very Large Telescope (VLT). The flux-calibrated, 
medium-resolution X-shooter spectrum (300 - 2500 nm) includes over 300 emission lines, but no 
(photospheric) absorption lines and is consistent with a reddened disk spectrum. Among the emission 
lines are three hydrogen series and helium lines, both permitted and forbidden metal lines, and CO 
first-overtone emission. A representative sample of lines with different morphologies is presented. 
The Ha and Ca II triplet lines are very strong, with profiles indicative of outflow and - possibly - 
infall, usually observed in accreting stars. These lines include a blue-shifted absorption component 
at ~ —125 km s _1 . The He I and metal-line profiles are double-peaked, with a likely origin in a 
circumstellar disk. The forbidden lines, associated with outflow, have a single blue-shifted emission 
component centered at —125 km s _1 , coinciding with the absorption components in Ha and Ca II. 
SINFONI H- and K-band integral-field spectroscopy of the cluster environment demonstrates that 
the [Fe n] emission is produced by a jet originating at the location of 08576nr292. Because the 
spectral type of the central object cannot be determined, its mass remains uncertain. We argue that 
08576nr292 is an intermediate-mass YSO with a high accretion rate (M acc ~ 10" 5 - 10" 6 M Q yr _1 ). 
These observations demonstrate the potential of X-shooter and SINFONI to study in great detail an 
accretion disk-jet system, rarely seen around the more massive YSOs. 

Subject headings: stars: formation — stars: mass-loss — stars: pre-main sequence — ISM: jets and 
outflows 



1. INTRODUCTION 



The formation process of massive stars is poorly un- 
derst ood (for a recent review, see Zinnecker & Yorke 
2007|. Besides complicated theoretical aspects, obser- 



vations of forming massive stars are limited due to their 
short formation timescale, the strong obscuration of their 
birth places, and their relatively small number. A cur- 
rent key question is whether massive stars form through 
accretion in a way similar to low- and intermediate-mass 



stars (e.g. 


Shu 


1977 


Palla & Stahler 


1993 


Sonnhalter 




Kuiper et al. 2010f, or that 



pressure and ste 



Yorke & 



lar wind prevent the formation of mas- 



sive stars, calling for alternative formation scena rios such 
as stellar mergers in dense and young c lusters ( |Bonnell 
et al.||1998| |Baumgardt fc Klessen||2010[ ). The detection 
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of the signatures of formation, i.e. circumstellar disks 
and collimated jets, around massive young stellar objects 
(MYSOs) would provide an important step in addressing 
this question. 

Up to now observational evidence for the presence of 
disks and/or collimated outflows has be en secured for 
abou t a dozen of (candidate) MYSOs (e.g. 
2007| ISandell fc Wright||2010t IKraus et a 
et al.[|2010| |Guzman et al.|[2010| |Zapata et"at]|2010p 



1'he current mass of these objects is often unclear, and 
their identification as a MYSO is mainly based on their 
luminosity. Most of them have a corresponding ZAMS 
mass less than 20 M (i.e. spectral type later than 08). 

We have initiated an observing campaign to extend the 
spectral coverage of a sample of MYSO candidates as far 
to the blue as possible, to better determine their pho- 
tospheric properties, to study the onset of stellar winds, 
and to characterize the physical structure of their ac- 
cretio n disks. A sample of MYSO candidates (Bik et al 



2006| hereafter B06) and young OB stars ( |Bik et al.|20D5t 
hereafter B05) has been obtained from a near-infrared 
survey of southern star-forming regions including ultra- 
compact Hn (UCHu) regions (Bik 2004). Among these 
is the region RCW 36, centered on the IRAS point source 
08576-4334, with the characteristic infrared colors of an 
UCHu. 



This point source i s associated with Vela Cloud C (Ya- 



maguchi et al. 1 999), for which th e estimated distance 



is 0.7 ± 0.2 kpc (|Liseau et al.|[l992[ ). The kinematic dis 
tance, however, from the rotation curve of Brand fc Blitz| 
( 1993 1 and a Local Standard of Rest velocity {vlsr) of 
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Obs.# 


date 


exp.time 


nod throw 


pos. angle 


I 


22-02-10 UT 5 h 06 


2x900 s 


2" 


19.52° (1) 


II 


22-02-10 UT 5 h 36 


2x900 s 


2" 


19.52° (1) 


III 


23-02-10 UT 4 h 50 


2x600 s 


5" 


50.26° (2) 



TABLE 1 

List of the observations carried out with X-shooter. The position 
angle is the slit angle on the sky from north to east. The numbers 
between parentheses refer to the slit orientations depicted in 



7.5 km s (Bronfman et al. 1996) is 2.2 kpc. In this 
work a distance of 0.7 kpc is adopted, consistent with 
the spectroscopic parallax obtained for OB stars in the 
cluster (B05). This region includes the MYSO candidate 
08576nr292. Based on J-K photometry, B06 estimate the 
spectral type of the central object at mid-B. 08576nr292 
has a strong infrared excess, and its K-band spectrum 
includes a broad Br7 emission line and CO ba nd-head 



emission, characteristic featu res of MYSOs (cf. Hanson 
et al.||1997| jBlum et ;il.||2()() 1|. The CO band-head emis 
sion can be modeled by a Keplerian rotating, circumstel- 



lar disk ( |Bik fe Thi|2004[ ). B06 prop ose that such a disk 
likely is a remnant accretion disk. [Wheelwright et ah] 
(2010) derive an inclination angle of i = 17.8^ U 4 S based 
on the CO-profiles, which means that the disk is seen 
almost pole-on. 

In the following we present the optical to near-infrared 
spectrum of 08576nr292 obtained with the new spectro- 
graph X-shooter, as well as near-infrared integral-field 
spectroscopy of its environment with SINFONI, both 
mounted on the ESO Very Large Telescope (VLT) at 
Paranal, Chile. The combination of these data sets re- 
veals important new information on the circumstellar en- 
vironment of this MYSO candidate. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. VLT '/ X-shooter 

The MYSO candidate 08576nr292 (RA (2000.0) 
08 h 59'21.7"; Dec (2000.0) -43°45'31.0") was observed as 
part of the X-shooter GTO program "Probing the ear- 
liest evolutionary phases of the most massive stars" in 
February 2010. A selection of the most important spec- 
tral features is presented in Fig. [T] Each observation 
comprised two separate exposures at different nodding 
positions on the slit (see Table [l]). In the NIR arm the 
exposure time was split into detector integration times 
(DITs) of 50s each. 

The X-shooter wavelength range is 300-2500 nm, split 
in three arms using dichroics: UVB (300 — 560 nm), 
VIS (550 - 1020 nm) and NIR (1000 - 2500 nm). For 
a detailed d escription of the instrument, see D'Odorico 



et al. (|2006j) . Given the strong reddening of the source, 
HnfeenFsIit widths were used in the UVB (1.0"), VIS 
(0.9") and NIR (0.4") arms. The resolving power was 
calculated from wavelength calibration frames and av- 
erages to 5100, 8800, and 11300 (corresponding to 59, 
35 and 27 km s _1 per resolution element) in the three 
arms, respectively. The V-band seeing varied between 
0.6" and 0.9" during the observations. A telluric stan- 
dard star (HD80055, A0V) was observed immediately af- 
ter 08576nr292. The spectrophotometric standard star 
GD71, a DA white dwarf, was observed in twilight. 



The ESQ X-shooter pipeline ve rsion 1.0.0 (Goldoni 
et al. 2006 Modigliani et al. 2010) was used to obtain 
reduced 2D spectra. From these 2D-spectra, ID-spectra 
were extracted and normalized to the continuum using 
our own routines. For the flux calibration we used a flux 
table for GD71 in the range 300 - 1000 nm, and scaled 
a Rayleigh- Jeans function to the 2MASS magnitudes for 
HD80055 in the range 1000 - 2500 nm. We estimated slit 
losses based on the seeing conditions. The derived J-, H- 
and K-band fluxes are within 30% of the 2MASS values: 
J = 11.82 ± 0.03, H = 10.38 ± 0.04, K = 9.32 ± 0.04. 

2.2. VLT/SINFONI 

Observations of the cluster environment of 08576nr292 
were made with the integral field spectrog raph SINFONI 
( |Eisenhauer et al.|2003{ Bonnet et al.|2004| . It produces 
an Ti- and K-band spectrum at every pixel (0.125" x 
0.125") in its field-of-view (8"x 8") at R - 1,500. The 
SINFONI data were collected on March 5, 7 and 19, 2007. 
Each observation was accompanied by the observation of 
an early B-type standard star observed immediately after 
the science observations, and used for flux calibration and 
telluric absorption correction. 

Different pointings in a raster pattern were used to 
cover a l arge area surro unding 08576nr292 (for more de- 
tails, see Bik et al. 2010 ). The DIT was 30 s per pointing. 
With every point being covered twice, this results in a 
total on-source integration time of 60 s. Sky frames were 
taken every 3 minutes using an empty area on the sky. 
The K-band seeing was 0.4" for the first two observing 
nights (part A in Fig. [2]). During the last night (part B), 
the seeing was 0.7". The south west corner of the image 
was not observed. 

The data were reduced using the SINFONI SPRED 
pipeline (ver sion 1.37) developed by the MPE SINFON I 



consortium (Schreiber et al. 



2004 Abuter et al. 2006) 



Velocity maps were created by fitting a gaussian profi 
to every spatial pixel in the data cubes. The velocities 
were converted to the LSR. 

3. RESULTS 
3.1. X-shooter spectrum of 08576nr292 

In the X-shooter spectra, the signal to noise ratio 
(SNR) is about 80 at 1200 nm and 50 at 650 nm, respec- 
tively, and less than 10 below 550 nm due to interstellar 
extinction. Despite the high SNR, no photospheric ab- 
sorption lines are detected. The spectrum contains over 
300 emission lines, almost all of which could be identi- 
fied. A representative sample of these lines highlighting 
the different morphologies is displayed in Fig. [T] The 
velocity scale is corrected for the Earth's motion with 
respect to the LSR. 

The hydrogen Balmer, Paschen and Brackett series are 
in emission, as well as many different metal lines. Close 
to half of all the identified lines are iron lines. Some 
variation in the peak-to-peak flux ratio was observed in 
some Fe I lines between observations I /II and III. 

(i) Double- and single-peaked emission lines: The 230 
permitted emission lines are predominantly from H I, 
He I, N I, O I, Mg I, Si I, Ni I, Fe I and Fe II. About 50 
of these lines (such as most He I and Fe I lines) are dou- 
ble peaked with a peak separation of 60 to 100 km s _1 
(Fig.[T^). The double-peaked signature is consistent with 
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Fig. 1. — A representative sample of the 300+ normalized spectral lines and their various profiles, detected by X-shooter in observation 
III of 08576nr292. (a) Double-peaked and single-peaked lines; (b) Emission lines with a blue-shifted absorption component at —125 km s 
and a 'red shoulder' (Ca II). The central absorption in Ha is due to nebular contamination; (c) Forbidden transitions blue-shifted at 
— 125 km s^ 1 ; (d) Peculiar profiles. The spectra are normalized. The rest wavelengths (nm) are in air. 



formation in a Keplerian rotating circumstellar disk, as 
suggested by the modelling of the CO band-he ads at 
2.3 /im present in the spectrum ( |Bik fe Thi||2004[ ). 

(ii) Emission lines with signatures of outflow and in- 
fall: There are a number of lines with strong, single 
peaks centered around the systemic velocity. Examples 
of these are the lower transitions in the Paschen and 
Brackett series, and the O I lines (Fig. [1^,). By far the 
strongest lines are Ha and the Ca II triplet (Fig. [T]d), 
which include a narrow absorption feature, blue-shifted 
to —125 ± 6 km s _1 , a signature of outflow. The Ca n 
lines have a red shoulder, a possible signature of infall 
onto the central star. A similar profile is seen in Pa/3. 
Note that the apparent double peak in Ha is an arte- 
fact of the subtraction of the nebular component, which 
varied over the slit. 

(Hi) Forbidden lines: A total of 67 identified forbid- 
den lines are detected, 46 of which areJFe n]; a repre- 
sentative sample is displayed in Figure [it. These lines 
have a single peak centered at a blue-shifted velocity of 
— 125 ± 8 km s _1 . This is the same velocity shift as 
that of the absorption features in Ha and Ca II. The 
peaks are slightly asymmetric with an extended wing to- 
wards the red and a sharp blue cutoff. The [O i] and 
[S n] lines have a very weak red-shifted counterpart at 
+140 ± 15 km s _1 . The [O i] doublet lines at 630 nm are 
the strongest of the forbidden transitions. These lines, 
as well as the [S n] lines, include a component at a low 
blueshifted velocity (~ —10 km s _1 ), that ma y be asso 



ciated with a hot wind close to the star (cf. Hartigan 



et al. 1995) 



b 'or bidden line emission is produced by collisionally ex- 
cited ions in a low-density medium. In this case, the two 
components at blue- and redshiftcd radial velocities sug- 
gest a bipolar outflow. The [Ca Ii] doublet at 729 / 732 
nm (Fig. [l^i) are the only detected forbidden lines that 



do not have a blue offset, but a double-peaked signature. 

(iv) Peculiar profiles: There are two peculiar line pro- 
files that do not fit in any of the above descriptions 
(Fig. fTH) . The first is the Na I A589 nm doublet, which 
includes an interstellar absorption feature, and a blue- 
shifted absorption component centered at —125 km s _1 . 
This absorption component is broader than that seen in 
Ca II and Ha. Another peculiar profile is seen in the He I 
A1083 nm line, which has a deep and narrow central ab- 
sorption component and a broad blue-shifted absorption 
trough extending from —200 up to —450 km s _1 . 

3.2. SINFONI data 

In Figure [2^,, a line map is plotted of the environment 
of 08576nr292 with tracers of ionized gas (Br7 A2166 nm 
in red), shocked or UV excited gas (H2 A2121 nm in 
green) and low-density gas ([Fe 11] A1644 nm in blue). 
A collimated bipolar outflow, from now on referred to 
as 'jet', is clearly detected in [Fe 11], originating at 
08576nr292. The jet is also detected in Br7 with a similar 
velocity structure. 

The velocity map (Fig. [2jj) shows that both lobes of 
the jet are moving away from the source. The velocity of 
the blue-shifted lobe corresponds to the peak velocity of 
the forbidden lines and the absorption features in Ha and 
Ca 11. In both lobes, some regions peak up to 200 km s _1 . 
Relatively high velocities are reached in those parts of the 
jet where the flux is relatively high (indicated by white 
contours). Careful inspection, however, reveals that the 
velocity peaks slightly closer to the source than does the 
flux. This suggests that fast-moving material is plough- 
ing into slower moving material. This is commonly seen 
in bo w-shocked jets that reveal an episod ic accretion his- 
tory ( |Raga et al.||1990| |Arce et al.||2007[ ). 

It is clear that the acceleration ot the jet material oc- 
curs close to the star (< 0.5" = 350 AU). The projected 
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Fig. 2.— (a) Detail of the SINFONI image of IRAS 08576-4334 (size: 57" X 43"). Color channels trace: H 2 (green), Br7 (red) and 
[Fe II] (blue). The results from the different observing nights, as described in section |2.2| are separated by the orange dotted line. The 
X-shooter slit orientations are depicted (slit size: 1" X f 1"). (b) Velocity offset map (size: 25" X 18") of the [Fe n] A1644 nm line. The 
color gradient corresponds to the offset of the peak with respect to the «lsr- Grey color represents areas where no [Fe II] emission was 
detected. The white contours represent the line flux, with logarithmic contour values between l.OxlO 4 and 2.3x10 s erg s — 1 cm -2 /im -1 
sr — 1 . The black cross indicates the position of 08576nr292, where [Fe II] emission was blended with Br-12; this area is also colored grey. 
We note the presence of another jet system in the NE corner of the field of view. 



lengths of the lobes are 9,000 AU (blue-shifted south-east 
lobe) and 3,000 AU (red-shifted north-west lobe). Prac- 
tically no H2 emission is detected throughout the jet, but 
this might be due to the short exposure time. 

The jet does not seem to end in a bow-shock. The 
end of the receding (NW) lobe may be obscured by fore- 
ground dust and at the end of the approaching (SE) lobe 
the [Fe n] emission changes velocity (the green-colored 
area in Fig. [2Jd). It is unclear whether this material is 
related to the jet, or part of another structure in the ISM. 

4. DISCUSSION 
4.1. Evidence for a disk-jet system 

The different line profiles described above and the ex- 
tended emission morphology show that 08576nr292 is a 
source with multiple components. A sketch of our in- 
terpretation of the system is depicted in Figure 3l The 
typical line profiles that we associate with the different 
components are displayed. 

The circumstellar disk is the probable origin of the 
many double-peaked emission lines. The double-peaked 
morphology likely depends on line strength (and thus the 
location and extent of the line- forming region), as indi- 
cated by the hydrogen spectrum, where only the higher 
transitions (3 — > 10 and up) show double-peaked emis- 
sion lines, although not very pronounced. The lower ion- 
ization species (Mg I, Fe i) typically have small peak sep- 
arations, consistent with being formed in the outer part 
of a Keplerian disk, shielded from the ionizing flux. He I 
and Fe II have a larger peak separation, thus presumably 
formed in a more rapidly rotating inner disk region. The 
strength of some of the single-peaked emission lines as 
opposed to the double-peaked lines could point out that 
they are formed over a more extended region of the Ke- 
plerian disk. Additional contributions to these lines may 
come from an accretion column and an ionized disk wind. 

The broad blue-shifted absorption in He 1 A 108 3 nm 
is an indicator of outflow (Kwan & Fischer 20101. We 
note that absorption over the velocity range up to 



—450 km s _1 is not seen in any other line. The pres- 
ence of a collimated bipolar jet is unambiguously estab- 
lished by the SINFONI observations. From this we infer 
that the forbidden emission lines with a radial velocity 
~ 125 km s _1 are formed in the jet. 

The P-Cygni like profiles of Ha, Ca 11 and Na 1 arc 
similar to those seen in some accreting Herbi g Ae/Be 



(HAeBe) and "class ical" T Tauri stars (CTTS, |Hamann 
& Persson 1992a|b[ ). Just as in 08576nr292, in some of 
these stars the Ca 11 triplet lines are very strong and their 
flux ratios indicate saturation. Furthermore, in CTTS, 
these lines are known tracers of t he accretion-infall zone 
(e.g. |Muzerolle et al.p998||205l] ). 

The coincidence ot the blue-shifted absorption features 
in these lines and the peaks of the forbidden emission 
lines suggests that they trace the same region. The ab- 
sorption may take place in the jet or in the disk wind, 
which would then be the launch region of the jet. In 



the CTTS models oflLima et al. (20101 a very similar 
Ha profile is produced in an accretion column in combi- 
nation with an ionized disk wind. 08576nr292 could be 
accreting in a similar way. 

4.2. Evolutionary status and mass 

The presence of a jet and the accretion morphology of 
the line profiles indicate that 08576nr292 is a YSO with 
an active accretion disk, contrary to the earlier interpre- 
tation that the circumstellar disk is a remnant passive 
disk (B06). Its mass and accretion rate, however, remain 
uncertain. The absence of photospheric lines prevents a 
spectral classification. The shape of the spectral energy 
distribution (SED) includes information about the mass 
and luminosity of the source. 

T he SED has to be cor rected for interstellar extinction 
(cf. Cardelli et al.|[l989| , which is a main source of un- 
certainty^ B06 determine Ay = 12 mag based on the 
spectral classification of three O stars in the same clus- 
ter, but the reddening is known to vary locally. This is 
underlined by the fact that the diffuse interstellar bands 
(DIBs) observed in these O stars (Ellerbroek et al., in 
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Fig. 3.— Sketch of the 08576nr292 disk-jet system. The typical 
line profiles associated with the different components are indicated. 
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Fig. 4. — SED as observed by X-shooter (black), de-reddened 
with Ay = 8 mag (orange). Furthermore plotted are a power law 
fitted to the dereddened SED (/3 = -0.6, green), |Kurucz| ||1993 ' 
models for AOV and A5V stars with ZAMS radii (2.5 and 1.7 Hq 
at d = 0.7 kpc (blue). Indicated with red squares are the 2MASS 
fluxes. 



prep.) are lacking in 08576nr292. However, the limited 
visibility of the receding jet lobe suggests that it is not a 
foreground object. A relatively low extinction could be 
explained by the fact that the jet has carved out a 'hole' 
in the ambient molecular cloud. 

With our current knowledge, we make a refined es- 
timate of Ay assuming that Ry — 3.1 and d = 
0.7 kpc. The absence of absorption features indicates 
that the photosphere is heavily veiled due to the accre- 
tion. Adopting a power law for the SED (XF\ oc X 13 ) 
we e xpect —4/3 < < in the range 500 - 1000 nm 



(cf. Hillenbrand et al.||1992[ ). Dereddening the spec- 
trum to^fialcTiTrnsTrat SED leads to Ay = 8 ± 1 mag. 
The assumption that disk and/or accretion luminosity 
dominates the SED allows to determine a limit on the 
spectral type. As can be seen from Fig. |1J this would 
exclude (PMS equivalents of) stars with spectral type 
earlier than AOV. 

Several remarks should be made. Firstly, the optical 
SED might not be flat or slightly declining, but increas- 
ing due to a cold gas disk. However, the many double- 
peaked emission lines, and the strong Ha, suggest that a 
hot and (partly) ionized disk is present. Secondly, a pos- 
sible underestimate of the distance would imply a larger 
u pper bound on the stell ar luminosity. 



Hartigan et al. 



(1995) derive empirical relations be- 
tween the line fluxes of optical forbidden lines and the 
mass-loss rate through the jet for CTTS. Applying these 
relations to the de-reddened spectrum of 08576nr292, we 
determine a jet mass-loss rate in the order of 10 -8 - 
10~ 7 Mq yr _1 . The same study suggests that the mass 
accretion rate is typically 10 2 times the mass-loss rate, 
implying an accretion rate of 10 -5 - 10 -6 Mq yr . Only 
a few intermedi ate-mass stars have been observed with 
such high rates ( Garcia Lopez et al.||2006| ). In this case 
one would expect the disk to be optically thick up to 
distances close to the stellar surfa ce (up to a few stellar 
radii, cf. Hillenbrand et al.||1992|). This would be con- 



sistent with the disk dominating the SED up to optical 
wavelengths. 

These considerations suggest that 08576nr292 is a 
HAeBe-like, Lada Class II source, undergoing active ac- 
cretion. Because of its brightness, relatively low extinc- 
tion and low inclination angle, 08576nr292 provides a 
unique opportunity to study an active disk-jet system in 
great detail. Further X-shooter observations and mod- 
elling of the jet are planned, aimed to further explore the 
dynamics of this system. 
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Appendix: 
X-shooter spectrum of 08576nr292 



Observed 23-02-2010 UT 4h50 
Not corrected for telluric absorption lines 
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